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In order to resolve the controversy about the low density region of the phase diagram of the 4 He 
monolayer on graphite, we have undertaken a path integral Monte Carlo study of the system. We 
provide direct evidence that the low density monolayer possesses solid clusters and a low density 
vapor as opposed to the most recent proposal that the system is in a superfluid phase. We further 
establish that the rounded heat capacity peaks observed at low densities are caused by melting of 
such solid clusters and are not associated with the suggested superfluid transition. 



PACS numbers 67.70. +n, 67.40 Kh 

Monolayer helium adsorbed on graphite has long 
proven to be a fascinating system for studying the growth 
and behavior of a quantum film above an atomically or- 
dered, uniform substrate and has been used to investi- 
gate a number of nearly two-dimensional (2D) phenom- 
ena • A very prominent feature of this layer is a com- 
mensurate solid phase in which one-third of the available 
substrate adsorption sites are occupied For den- 

sities above the commensurate, the system is known to 
pass through a region of domain wall phases before form- 
ing an incommensurate triangular solid phase. On the 
other hand, the nature of the phase diagram below the 
commensurate density at low temperatures is less well es- 
tablished, with two competing pictures. One possibility 
is that the phase is a solid with vacancies || . At low tem- 
peratures, the vacancies coalesce, producing coexistence 
between a solid cluster and a vapor. More recently, it has 
been suggested that the solid melts if the density is de- 
creased, and there is a low temperature liquid phase 0. 
The first layer would then be a candidate in the search 
for a monolayer superfluid p|JlO|] . 

The solid cluster picture has been discussed by Ecke et 
al Q . They note that since the commensurate solid phase 
is in the same universality class as the three-state Potts 
model [O-O], then at lower densities the film should 
consist of a commensurate solid with vacancies. If the 
temperature is raised, the solid melts continuously. Low- 
ering the temperature causes the vacancies to coalesce 
(phase separate), a first order transition. The difference 
between the temperatures of these two transitions be- 
comes smaller as the density is lowered until they meet 
at a tricritical point. This point was determined to be 
at about 0.039 atom/A 2 and 1.3 K. Thus, in this picture, 
the monolayer consists of solid clusters surrounded by a 
low density vapor at low temperatures and densities. 

More recent experiments have questioned this conclu- 
sion. Grey wall and Busch || point out that the heat ca- 
pacity is not linear in density for the entire region below 
the commensurate density, as it must be for solid-vapor 
coexistence jjj. They instead propose that the system 
has a self-bound liquid phase at about 0.04 atom/A 2 . 



This conclusion is supported by 2D variational calcula- 
tions for helium Q that take substrate corrugations into 
account. The possibility of a first layer liquid is intrigu- 
ing since at low temperatures it would be a superfluid 
above a bare substrate, with no underlying "dead" layer 
of helium. However, direct measurements fl5| on the 
first layer detect no superfluidity. This negative result 
has been attributed to poor substrate connectivity. 

In order to resolve this controversy, we have under- 
taken a path integral Monte Carlo (PIMC) study of the 
low density first layer. This is the first attempt to di- 
rectly investigate this region by an exact, first principles 
method that treats the full quantum many-body prob- 
lem. From our calculations, we provide the first direct 
evidence that the low density monolayer consists of solid 
clusters. No liquid phase occurs, and so there is no pos- 
sibility for first layer superfluidity. We further establish 
that the rounded heat capacity peaks observed at low 
densities are caused by the melting of solid clusters and 
are not associated with a superfluid film, as has been 
suggested [||. 

Our PIMC calculations use realistic helium- helium |l6| ] 
and helium-substrate interactions. In order to include 
the effects of substrate corrugations, we use the full, 
anisotropic helium-graphite potential of Ref. |L7|]. For 
a general discussion of the PIMC method and its appli- 
cation to films, see Refs. [[l8],[l9|. In test runs, we deter- 
mined that an inverse temperature slice of r = 1 /200K~ 1 
was required to reach the desired accuracy using the semi- 
classical approximation for the required high tempera- 
ture density matrix at 200 K. The semiclassical approxi- 
mation allows substrate corrugations to be easily imple- 
mented. From the same test runs, we determined that an 
/ = 3 multilevel bisection was required. All of our calcu- 
lations are performed in a simulation cell with periodic 
boundary conditions and dimensions 25.560A x 22.136A 
that exactly accommodates the commensurate solid The 
number of particles ranged from 20 to 40, with 36 cor- 
responding to the commensurate density, p c — 0.0636 
atom/ A 2 . We also allowed for the possibility of particle 
permutations at low densities, but did not observe any. 
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It is essential that corrugations be included in the cal- 
culations because the commensurate solid phase is en- 
tirely the result of the substrate corrugations. A recent 
simulation for the helium monolayer ]20[ ] using the later- 
ally averaged graphite potential (nj finds that the equi- 
librium phase is a liquid. Full solidification does not oc- 
cur until the coverage is well above p c . At p c , the film on 
the featureless substrate is a compressed uniform liquid 
that is near the beginning of solid-liquid coexistence. 

Before presenting evidence that the first layer has solid 
clusters, we first wish to demonstrate that our simu- 
lation method can reproduce the commensurate solid 
phase, and that this phase exhibits melting-like behav- 
ior in agreement with experiment. We then investigate 
the low temperature phase diagram using the Maxwell 
construction. 
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FIG. 1. Distribution plots at the commensurate 
0.0636 atom/A 2 , for T=2.99 K (left) and T=4.0 K 
Filled circles indicate graphite adsorption sites. 
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density, 
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Figure [l] illustrates the V3 x \/3 commensurate solid 
phase and its melt using probability density contour 
plots. Raising the temperature from 3 K to 4 K causes 
melting, so that each adsorption site will, after a suffi- 
ciently long simulation run, have an equal probability of 
being occupied. Further evidence for solidification comes 
from static structure factors. Figure ^ shows our calcula- 
tions at and immediately below the commensurate solid 
density for the (01) scattering direction. The peaks at 
1.70 and 3.40 A -1 are the wave vectors expected for the 
first two Bragg scattering peaks for the -s/3 x V3 solid. 

The temperature dependence of the static structure 
peaks can be used to determine melting temperatures. 
Figure || shows S(k — 1.7 A~ 1 )/N for several densities. 
Melting is signaled by a drop in the average peak height 
and large statistical fluctuations in peak values. This 
is first observed at 2 K, 2.5 K, 3 K, and 3.33 K for 
0.0424, 0.0530, 0.0566, and 0.0636 atom/A" 2 , respec- 
tively. This density dependence of melting is in agree- 
ment with the experimental phase diagram, although our 
melting temperatures are somewhat higher than the ex- 
perimental values. Heat capacity measurements indicate 
the commensurate solid melts at 3 K, and the low density 
(< 0.045 atom/A -2 ) melting peaks are at about 1.5 K. 

We have also calculated the temperature dependence of 
the energy per particle for several densities. These pos- 
sess inflection points that produce heat capacity peaks 
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FIG. 2. Static structure factor S(k) for 0.0530 (diamonds), 
0.0566 (triangles) and 0.0636 atom/A 2 (squares) at 2.0 K. 
The particle numbers are are 30, 32, and 36, respectively. 



when the values are differentiated, signaling melting. 
Two sample calculations of the specific heat are shown 
in Fig. |^. At the commensurate density, we observe 
that melting occurs at about 3.5 K, somewhat above the 
experimental value but consistent with the static struc- 
ture calculations. As the density decreases, the melting 
temperature and peak height also decrease. At 0.0353 
atom/A 2 , the melting temperature is about 1.5 K, and 
the peak is much smaller and more rounded. 

The binding energy Eg of a single particle on the 
substrate may be easily calculated. We find Eg = 
— 145.48±0.21, which is comparable to, but slightly lower 
than, the estimated values of — 141. 75± 1.50 K from scat- 
tering 21 L and —142.33 ± 1.97 K from thermodynamic 



analysis |22[| . Subtracting Eb from the energy per parti- 
cle at p c , we find the solid binding energy is —2.21 ± 0.20 
K. This is lower than either the energy per particle of the 
2D liquid or solid determined in Ref. Q. 

Having demonstrated that our simulation method can 
reproduce known features of the monolayer, we now turn 
to the low density phase. This region is investigated by 
applying the Maxwell common tangent construction to 
the low temperature values of the total energy to iden- 
tify unstable regions at effectively zero temperature. The 
application of this method to a system with a constant 
volume and varying particle number is described in Refs. 
p9[ . The free energy at nonzero temperatures is not di- 
rectly accessible from PIMC calculations, but we may 
determine effectively zero-temperature energy values by 
a limiting process p9| . The total free energy and total 
energy are the same at zero temperatures. 

The coexistence region between two stable phases is 
characterized by an unphysical upward curvature of the 
total (free) energy. In the thermodynamic limit, the en- 
ergy values lie on a coexistence line. Such an unstable 



2 



0.50 



0.40 - 



0.10 - 




0.00 1 1 1 1 1 1 1 1 1 1 1 1 

0.8 1.2 1.8 2.2 2.8 3.2 

Temperature (K) 

FIG. 3. Temperature dependence of the peak height of 
the static structure factor. The coverages are 0.0424 (cir- 
cles), 0.0530 (filled diamonds), 0.0566 (triangles), and 0.0636 
atom/ A 2 (filled squares). 

region may be identified between zero coverage and p c 
in the total energy values shown in Fig. 0. We have 
verified that all these energy values have approached the 
zero-temperature limit within error bars. All interme- 
diate energy values are above the coexistence line. In 
the unstable region, the system can phase separate into 
a zero density vapor and a commensurates solid cluster. 
The energy increase results from the finite cost of creat- 
ing the phase boundary. 
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FIG. 4. The heat capacity at 0.0353 (filled circles) and 
0.0636 atom/A 2 (squares). The dashed line is a guide to the 
eye. The solid line is the measured specific heat at 0.0367 
taken from Ref. |j] 

Contour plots of the probability density, shown in Fig. 
^, provide direct evidence that both solid-vapor coexis- 



tence and solid phases with vacancies occur. At the lower 
temperature (left-hand side of Fig. ||), the vacancies have 
coalesced. Note that there is only one bubble in the left- 
hand side of Fig. |^ because of periodic boundary condi- 
tions. The holes move very slowly at this temperature, 
producing long equilibration times for condensation. We 
thus calculated the energy for this system with the vacan- 
cies initially separated and then initially condensed. The 
energy for condensed vacancies was lower. At higher tem- 
peratures, the vacancies acquire enough kinetic energy to 
leave the phase separated state and diffuse into the solid. 
As a result, vacancies can become isolated. This is illus- 
trated at 2.5 K in the right-hand plot of Fig. [| A series 
of probability distribution plots reveals that these vacan- 
cies move in the simulation, so the equilibration problem 
encountered at 1.0 K is not present at this temperature. 
We note that in our simulation there is still evidence of 
phase separation in contour plots at 2.0 K for the den- 
sity shown in Fig. [| while experimental results seem to 
indicate a transition at 1.5 K. We have plotted probabil- 
ity contours for densities as low as 0.0207 atom/A 2 and 
observe solid clusters at all densities. 

We have also attempted to place a vacancy in a solid 
cluster surrounded by vapor at 0.0424 atom/A 2 and 1.0 
K. The vacancy was spontaneously expelled from the 
cluster during thermalization, from which we conclude 
that the solid clusters cannot also contain isolated va- 
cancies when in equilibrium with the low density vapor. 
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FIG. 5. Total energy versus coverage. For clarity, the en- 
ergy values have been shifted by the line Neo, where N is the 
number of particles, and eo = —147.78 ± 0.08, the minimum 
energy per particle, eo occurs at the commensurate density. 
The dashed line is the gas-solid coexistence line. 

Finally, we wish to discuss the arguments of Greywall 
and Busch (GB) against solid clusters and in favor of the 
superfluid phase. Their primary objection to solid-vapor 
coexistence is that this should be signaled by linear heat 
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capacity isotherms for the entire region from zero cov- 
erage up to the commensurate density. Their published 
data shows that for temperatures from 0.2 K to 0.5 K, 
the isotherms are linear only between 0.025 and 0.060 
atom/A 2 . At 0.1 K, the upper endpoint is about 0.055 
atom/ A 2 . As a possible explanation, we suggest that the 
departure from linearity below 0.025 atom/A 2 is caused 
by the presence of multiple finite-sized clusters. At low 
densities, solid clusters nucleate around surface defects. 
Initially, there are many small metastable clusters with 
large perimeter-to-area ratios. Increasing the density in- 
creases the size of the clusters until the surface is covered 
by a few large solid clusters with negligible boundary ef- 
fects. Thus, the heat capacity exhibits linear behavior 
only after the solid clusters are sufficiently large so that 
the pcrimctcr-to-arca ratio is small. This presumably oc- 
curs for coverages above 0.025 atom/A 2 . GB have used 
a similar explanation in their arguments for solid-liquid 
and liquid-gas coexistence in regions that do not have 
linear isotherms. 



performed using the facilities of the Supercomputer Com- 
putations Research Institute at Florida State University. 
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FIG. 6. Probability distributions for 0.0566 atom/A 2 at 
T=1.0 (left) and 2.5 K (right). The filled circles give the 
locations of graphite potential minima. 

GB are lead to identify coverages near 0.04 atom/A 2 
as liquid based partly on simulation results for 2D he- 
lium on a flat substrate, which were the most relevant 
calculations then available. As GB note, the large peak 
associated with the melting of the uniform commensurate 
solid phase first emerges at 0.04 atom/ A 2 . 2D helium is 
a liquid near this density |23] , suggesting that first layer 
coverages below 0.04 may be liquid. Unlike the purely 
2D simulations, our calculations take the role of sub- 
strate effects into account. As we have shown, surface 
corrugations push the density of the energy minimum up 
from about 0.04 on a flat substrate to 0.0636 atom/ A 2 
and produce solidification. GB also show that their heat 
capacity results are in general agreement with a PIMC 
calculation for 2D superfluid helium , suggesting that 
there might be a superfluid transition in the first layer. 
We have shown in Fig. ^ that the rounded heat capaci- 
ties seen for low first layer densities are produced by the 
melting of a solid cluster and are not associated with a 
superfluid transition. 

This work was supported in part by the National Aero- 
nautics and Space Administration under grant number 
NAG3-1841. Some of the calculations for this work were 
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